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• Abstract 

• The interlace polynomials introduced by Arratia, BoUobas and Sorkin 

extend to invariants of graphs with vertex weights, and these weighted 

■ interlace polynomials have several novel properties. One novel property 

U' 

^ '. qiG) = q{G -a) + q{G^'' - b) + {{x - if - l)q{G^'' - a - 6) 

» . 

that lacks the last term. It follows that interlace polynomial computations 
can be represented by binary trees rather than mixed binary-ternary trees. 
Binary computation trees provide a description of q{G) that is analogous 
to the activities description of the Tutte polynomial. If G is a tree or 
forest then these "algorithmic activities" are associated with a certain 
kind of independent set in G. Three other novel properties are weighted 
' pendant-twin reductions, which involve removing certain kinds of vertices 

QQ , from a graph and adjusting the weights of the remaining vertices in such 

a way that the interlace polynomials are unchanged. These reductions 
allow for smaller computation trees as they eliminate some branches. If a 
graph can be completely analyzed using pendant-twin reductions then its 
interlace polynomial can be calculated in polynomial time. An intuitively 
\ pleasing property is that graphs which can be constructed through graph 

substitutions have vertex-weighted interlace polynomials which can be 
obtained through algebraic substitutions. 
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1 Introduction 

Motivated by problems arising from DNA sequencing and the properties of 
circle graphs of 2-in, 2-out digraphs, Arratia, BoUobas and Sorkin introduced 
a new family of graph invariants, the interlace polynomials, in [3llll[S]. These 
invariants may be defined either through recursive elimination of vertices or as 
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sums indexed by subsets of vertices [TJ [S] , much as the Tutte polynomial may 
be defined either through recursive elimination of edges or as a sum indexed by 
subsets of edges. EUis-Monaghan and Sarmiento |2D] have shown (among other 
results) that one of the one- variable interlace polynomials, qn, can be computed 
in polynomial time for bipartite distance hereditary graphs. Their proof depends 
on the corresponding result for the Tutte polynomials of series-parallel graphs, 
first proved in [26], and the fact that bipartite distance hereditary graphs are 
circle graphs of Euler circuits in medial graphs of series-parallel graphs. 

In this paper we discuss several useful features of interlace polynomials that 
have been modified to incorporate vertex weights. After defining the weighted 
interlace polynomials in Section 2, we observe that a simple adjustment of 
weights makes it unnecessary to have the third term in the fundamental recur- 
sion q{G) = q{G - a) +q{G''^ -b) + {{x-lf- l)g(G"^ - a - 6) of [S]. In Section 
3 we present reduction formulas that can be used to eliminate a vertex that is 
a twin of another, or pendant on another. These pendant-twin reductions are 
analogous to the series-parallel reductions of electrical circuit theory, in which 
two resistors wired in series (resp. parallel) are replaced by one resistor with 
R= Ri+ R2 (resp. = R^^ + R^^)- The pendant-twin reductions are used 
to extend the result of EUis-Monaghan and Sarmiento mentioned above to the 
two-variable interlace polynomials of looped, non-bipartite distance hereditary 
graphs. In Section 4 we show that \iG = H^K \sa. looped graph obtained using 
the composition construction of Cunningham [19] then q{G) is equal to the in- 
terlace polynomial of a suitably re- weighted version of K\ this generalizes results 
of [4] that describe the q^ polynomials of simple graphs constructed through 
substitution. Composition has proven useful in the study of circle graphs (see 
for instance [TTl [TH [T^), so it is not surprising to see it appear in the theory 
of the interlace polynomials. In Section 5 we discuss some elementary proper- 
ties of the unweighted q^ polynomial, focusing on simple (unlooped) graphs. 
In Section 6 we sketch a combinatorial description of the interlace polynomials 
of trees and forests introduced by Anderson, Cutler, Radcliffe and the present 
author in [5]. This combinatorial description bears a striking resemblance to 
the activities description of the Tutte polynomial, and in Section 7 we extend it 
to arbitrary graphs using "activities" defined with respect to recursive interlace 
polynomial calculations. We do not know whether or not these activities have 
a convenient combinatorial description in general. 

We should observe that the idea of using vertex weights for interlace poly- 
nomials has appeared before, though our implementation of the idea is different 
from those we have seen elsewhere. In |16) , Courcelle introduced a multivariate 
interlace polynomial that is more complicated than the polynomials we consider 
here, and involves assigning indeterminates to the vertices of a graph. He used 
monadic second-order logic to show that it is possible to compute bounded por- 
tions of this polynomial (and the entire unweighted interlace polynomial q) in 
polynomial time for graphs of bounded clique-width. This technique is quite 
general but involves large built-in constants, so for pendant-twin reductions 
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and compositions the formulas presented here are considerably more practical. 
Also, Blaser and Hoffmann [7j use the idea of assigning indeterminates to ver- 
tices, along with the adjunction of two types of pendant-twin vertices, to show 
that evaluating interlace polynomials is generally #P-hard for almost all values 
of the variables. 

In [3] Arratia, BoUobas and Sorkin observe that there is a natural (so natural 
it is "practically a tautology") correspondence between the Kauffman bracket 
of an alternating link diagram and an interlace polynomial of an associated 2- 
in, 2-out digraph. The situation is clear enough that we do not discuss it in 
detail, but it is worth mentioning that this correspondence may be extended 
to arbitrary link diagrams using vertex weights. The well-known relationship 
between the Jones and Tutte polynomials is similar, in that an edge-weighted 
or -signed version of the Tutte polynomial conveniently incorporates crossing 
information when dealing with non-alternating links [231 1251 127j . 

2 Expansions and recursions 

We recall terminology and notation of [5j. Graphs may have loops but not 
multiple edges or multiple loops. The rank r{G) and nullity n{G) of a graph 
G are those of its adjacency matrix considered over GF{2). If S" C V{G) then 
G[S] is the subgraph of G induced by S. In addition, we say a graph is {vertex-) 
weighted by functions a and /3 mapping V{G) into some commutative ring with 
unity R. For ease of notation we prefer to denote a weighted graph G rather than 
using the triple (G, ol, (3); even when two graphs differ only in their weights we 
will denote them G and G' rather than (G, a, (3) and (G, a' ,(3'). Also, if a graph 
is modified then unless otherwise stated, we presume that the weight functions 
a and (3 are modified in the most natural way. For instance, if a G V{G) then 
the other vertices of G have the same vertex weights in G — a as they have in 
G. If G is a weighted graph then the unweighted version of G is denoted G''; it 
has the same underlying graph and the trivial weights a = /3 = 1 G Z. 

Definition 1 If G is a vertex-weighted graph then the weighted interlace poly- 
nomial of G is 

SCV(G) sGS v^S 

Definition 2 If G is a vertex-weighted graph then the weighted vertex-nullity 
interlace polynomial of G is 

<in{g)= e (n"(^))(n/5("))(y-i)"^''^''"- 

SCV(G) seS v^S 

Definition 3 If G is a vertex-weighted graph then the weighted vertex-rank 
interlace polynomial of G is 
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The original, unweighted interlace polynomials of G are recovered by using 

G". 

The three weighted interlace polynomials of a particular weighted graph G 
may be substantially different from each other. Considered as functions defined 
on the class of all weighted graphs, however, the three polynomials are essentially 
equivalent: if G is obtained from G by re-weighting V{G) = {vi, . ..,«„} using 
the indeterminates in the polynomial ring Z[ai, .., an, Pi, (3n], then the value 
of any one weighted interlace polynomial on G determines the values of all 
three polynomials on all weighted versions of G. For instance, q{G) can be 
obtained from qN{G) by substituting (x — 1) • a{vi) for a^, /3{vi) for Pi and 
1 + for y. (In contrast, the functions defined by the three unweighted 
interlace polynomials are rather different from one another, as there are many 
pairs of graphs distinguished by g^j but not by [5].) We mention all three 
weighted polynomials simply because one of them may be more convenient for 
some purposes than the others. Most of our results are stated for q because it 
specializes to the others most readily. 

In fact any one of q{G), q]^{G), (?_r(G) determines G up to isomorphism, 
because the vertex weights identify the contribution of each S C V{G). The 
looped vertices of G appear in the 1-element subsets S of rank 1 (nullity 0), 
some pairs of adjacent vertices of G appear in the 2-element subsets S that 
contain at least one unlooped vertex and have rank 2 (nullity 0), and the other 
pairs of adjacent vertices of G appear in the 2-element subsets S that contain 
two looped vertices and have rank 1 (nullity 1). Moreover the same comment 
applies even if only one of a, /3 is nontrivial, i.e., if = 1 or a = 1. This might 
make it seem possible to simplify our discussion by considering only a or only 
P, but we prefer to use both weights because they produce especially easy-to- 
read formulas, as noted in the discussions of Theorems 171 and fT2l below. Setting 
one weight identically equal to 1 would not simply leave us with the other 
weight; in essence, the remaining one would replace the ratio of the original 
two. Consequently, using only a or only P would entail unnecessary algebraic 
complications and losses of generality, because many formulas would require 
division. For instance, with P = 1 Proposition [4] would state that replacing 
a{a) with a' (a) = ria{a) / {ri + results in a graph G' with (ri + r2)q{G') = 
TiqiG) + r2q{G — a). The special case ri = — r2 would require a separate 
statement, and of course /3(a) =0 would be ruled out. 

Our first proposition follows immediately from Definition [TJ 

Proposition 4 Suppose a G V{G) and ri,r2 G R- Let G' he obtained from G 
by changing the weights of a to a' [a) = ria(a) and P'{a) — riP{a) + r2. Then 
q{G') = r^q{G)+r2q{G-a). 

A fundamental property of the interlace polynomials is that they can be 
calculated recursively with the local complementation and pivot operations used 
by Kotzig [24|, Bouchet OHi] and Arratia, BoUobas and Sorkin OHUl]. 
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Definition 5 (Local Complementation) If a is a vertex ofG then G° is obtained 
from G by toggling adjacencies {x,y} involving neighbors of a that are distinct 
from a. 

Definition 6 (Pivot) If a and b are distinct vertices of G then the graph G""^ 
is obtained from G by toggling adjacencies {x,y} such that x,y (f. {a,b}, x is 
adjacent to a in G, y is adjacent to b in G, and either x is not adjacent to b or 
y is not adjacent to a. 

Note that local complementation includes loop-toggling (when x = y is & 
neighbor of a distinct from a) but pivoting does not. 

Theorem 7 If G is a weighted graph then q{G) can be calculated recursively 
using the following properties. 

(a) If a is a looped vertex then 

q{G) = (3{a)q{G - a) + a{a){x - l)q{G'' - a). 

(b) If a and b are loopless neighbors then 

q{G) = (3{a)q{G-a)+(3{b)q{G''^-b) + {a{a)a{b){x-lf-l3{a)(3{b))q{G''^-a-b). 

(c) If G has no non-loop edges then q{G) is 

n {aiv)iy - 1) + Piv))]- I n («(«)(-x-l)+/3M) 

^unloopod veV(G) ) yioopcd v&(G) 

Proof. The proofs of parts (a) and (b) of Theorem [7] are essentially the same 
as the proofs of the corresponding formulas for the unweighted two-variable 
interlace polynomial [5] . Indeed, if we read a as "includes" and /3 as "excludes" 
then the weighted formulas serve as mnemonic devices to recall the proofs. 

For instance part (a) is proven as follows. The term f3{a)q{G — a) reflects the 
fact that if a ^ 5 C V{G) then G[S] (G - a)[S], so the contributions of S to 
q{G) and q{G — a) differ only by a factor (3{a). The term a{a){x — l)q{G"' — a) 
reflects the fact that if a G S* C V{G) then 

r{G[S]) - r I 1 Mn M,2 | - 1 + r (^^^JJ = 1 + r{{G^ - a) [5 - a]), 

where bold numerals represent rows and columns and Aff^ differs from Mn 
in every entry. Consequently the contributions of S to q{G) and S — {a} to 
q{G'^ — a) differ by a factor a{a){x — 1). 

The more complicated formula of part (b) reflects, among other things, the 
fact that each subset S C V(G) with a^b ^ S contributes to all three terms; the 
last two contributions cancel each other. 

Part (c) follows directly from Definition [TJ It is not technically necessary to 
mention graphs with loops in (c), as loops can always be removed using part 
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(a). However applying part (a) to completely disconnected graphs is obviously 
inefficient. ■ 

Other properties of the unweighted interlace polynomials also extend natu- 
rally to the weighted polynomials. For instance, the second and third parts of 
Theorem [8] below extend two properties discussed in Section 3 of [5] . The first 
part will be useful in Section 4, where we discuss substituted graphs. 



Theorem 8 (a) If a E V{G) is loopless then 

q{G) - (3{a)q{G -a) = q{G-) - f3{a)q{G'^ - a). 

(b) If a,b € V{G) are loopless neighbors then 

q{G - a) - (3{a)q{G - a - b) = q{G'''' -a)- /3(a)g(G"''' -a-b). 

(c) IfG is the union of disjoint subgraphs Gi and G2 then q{G) — q{Gi)q{G2)- 

Proof. To prove (a), note that if a e 5 C V{G) then the adjacency matrix 
of G[S] may be represented by 

^0 

Adding the first row to each of those in the second group results in 

^0 
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the adjacency matrix of ^[S]. Definition [T] then tells us that 
g(G) - P{a)q{G - a) 

aeS<ZV(G) seS v^S 

aescy(G) ses v^s 

= q{G-)-f3{a)q{G- -a). 

The proofs of (b) and (c) are essentially the same as the proofs of the corre- 
sponding results in [5J. ■ 
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Many properties of the unweighted interlace polynomials extend directly to 
the weighted polynomials, as we see in Theorems [7] and [8] It may be a surprise 
that some properties of the unweighted interlace polynomials can be signifi- 
cantly simplified using vertex weights. For instance, consider a computation 
tree representing a recursive implementation of Theorem [71 The tree has two 
branches for each application of part (a), three branches for each application 
of part (b), and a leaf for each application of part (c). As noted in the proof 
of Theorem [71 the three terms of part (b) all incorporate contributions from 
the same subsets S C V{G). The recursive computation calculates these same 
contributions three times, on separate branches. This inefficiency can be elim- 
inated by rephrasing part (b) of Theorem [7| so that no three-fold branches are 
necessary. 

Corollary 9 If a and b are loopless neighbors in G then q{G) ~ j3{a)q{G — a) -\- 
a{a)q{{G°'^ — b)'), where (G"'' — 6)' is obtained from G°'^ — b by changing the 
weights of a to a' {a) = f3{b) and 13' {a) = a{b){x — 1)^. 

Proof. Let (G"'' — 6)" be the graph obtained from G""^ — & by changing the 
weights of a as in Proposition [H with ri ~ P{b) and r2 — a{a)a{b){x — 1)^ - 
/3(a)/3(6). Then a" (a) = a(a)/3(6) and/3"(a) = a{a)a{b){x-l)^ , so g((G'^''-6)") 
= a{a)q{{G°''' - b)'). On the other hand. Proposition H tells us that 

qiiG"^ - b)") = f3ib)q{G'"' -b) + {a{a)aib){x - if ~ /3(a)/3(6))g(G'^'' -b-a). 



If a'{a) — (3{b)/{x — 1) and f3'{a) — a{b){x — 1) are used instead of the 
weights given in CoroUarylH then the resulting formula q{G) = f3{a)q{G — a) + 
a{a){x — l)q{{G°''' — b)') still has only two branches, and bears an interesting 
resemblance to part (a) of Theorem [71 However using division to define a' (a) 
may occasionally cause some algebraic difficulties, e.g., it complicates the eval- 
uation at 2: = 1, and it prohibits the use of rings in which a; — 1 is a divisor of 
zero. 

Corollary [9l is one of several results in which vertex weights allow us to 
extend properties of the unweighted version of ^at, seemingly the simplest kind 
of interlace polynomial, to the other interlace polynomials. In this instance the 
extended property is that of possessing a recursive description represented by 
a binary computation tree. Another result of this type is Corollary [TOl which 
extends Remark 18 of [4]: if a and b are loopless neighbors in G then the 
unweighted vertex-nullity interlace polynomials of G and G"** are the same. 

Corollary 10 Let a and b be unlooped neighbors in G, and let [G""^)' be the 
weighted graph obtained from G""^ by changing the weights of a and b to a' (a) = 
I3{b), I3'{a) = a{b){x - if, a'{b) = l3{a) and l3'{b) = a{a){x - if. Then 
ix-lfq{G) = qiiG->'y). 

Proof. Applying Corollary [9l to (G°^)', with the roles of a and b reversed, 
tells us that g((G'^'')') = /?'(fe)g((G'^'')' - 6) + a'(6)g((((G"'')')"'' - a)'). Observe 
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that (((G"^)')"'' - a)' has the underlying graph (G"'')"^ - a = G - a, and differs 
from G — a only in the weights of 6, which are given by a"{b) — f3'{a) — a{b){x — 
1)2 and/?"(6) =a'(a)(a;-l)2 = P{b){x-lf. Consequently a' (6)g((((G°^)')"'' - 
a)') = a'{h){x - \fq{G - a) = /3(a)(a; - \fq{G - a). As 13' {b)q{{G''^)' -h) = 
a{a){x — 1)2(7((G°^ — 6)'), the result follows directly from Corollary^ ■ 

If we are confident that division by x — 1 will cause no trouble, Corollary 
[To] may be restated with a simpler conclusion: using a' {a) = f3{h)/{x — 1), 
(3'{a) = a{b){x - 1), a'{b) = (3{a)/{x - 1) and = a{a){x - 1) yields 

g(G) = g((G-^)'). 

3 Pendant-twin reductions 

Other novel properties of the weighted interlace polynomials arise from a general 
observation. Suppose a, 6 e V{G) happen to have the property that for 5* C 
y(G), the rank and nullity of S are determined by the rank and nullity of 
S — {b}, perhaps in different ways according to which of a, b is contained in S. 
Then it may be possible to adjust a{a) and /3(a) in such a way that the weighted 
interlace polynomial of the weight-adjusted version of G — 6 incorporates all the 
information in q{G). The simplest instance of this observation occurs when a 
and b give rise to identical rows and columns in the adjacency matrix of G. 

Definition 11 Two vertices a,b of G are identical twins if (i) either they are 
looped and adjacent or they are unlooped and not adjacent, and (ii) they have 
the same neighbors outside {a, 6}. 

Identical twins are called clones in [7], and unlooped identical twins are 
called duplicates in [4J and false twins in |,20j . We prefer the present terminology 
because the adjective false seems inappropriate, and because we do not know 
what we would call non-identical duplicates or clones in Definition 7 below. 
The following result extends Proposition 40 of 0], Proposition 4.14 of [20] and 
Section 3.1 of 7j to vertex- weighted graphs. 

Theorem 12 Suppose a and b are identical twins in G. Let G' be the graph 
obtained from G~b by changing the weights of a: P'{a) — l3{a)(3{b) and a'(a) = 
a{a)(3{b) + a{a)a{b){y - 1) + /3(a)a(&). Then q{G) = q{G'). 

Proof. If a ^ S' C V{G') then r{G'[S]) = r(G[5]), because the adjacency 
matrices are the same. If a e S' C V{G') then r(G'[5]) = r(G[5]) = r(G[5U{6}]) 
= r{G[{S — {a}) U {6}]), because the only difference among the adjacency ma- 
trices is that the matrix of G[S U {6}] has two identical rows and columns, 
corresponding to a and b. As mentioned above, it is helpful to read a as "in- 
cludes" and /3 as "excludes" so that (for instance) the appearance oi (3{a)a{b) in 
a' {a) indicates that if a G 5 C V{G') then the summand oiq{G') corresponding 
to S includes the summand of q{G) corresponding to [S — {a}) U {6}. ■ 

Theorem 1121 has the following inductive generalization. Suppose fc > 1 and 
a — bojbi, ...,bk are identical twins in G. If G' is the graph obtained from 
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G-b 



•1 



bk by changing the weights of a to 



a (a) 



k 



and /3'(a) 



i=0 



then q{G) = q{G'). We refer to the process of combining several identical twins 
into a single vertex as an identical twin reduction no matter how many identical 
twins are combined. 

Theorem 13 Suppose b is an unlooped degree-one vertex pendant on a. Let G' 
be the graph obtained from G^b by changing the weights of a: a' (a) = a{a)j3(b) 



and 13' {a) = a{a)a{b){x - If + (3{a)a{b){v - I) + I3{a)li{b). Then q{G) =q{G'). 

Proof. Ua e S C V{G') then r{G'[S]) = r(G[5]), and if a ^ 5 C V{G') 
then r{G'[S]) = r{G[S U {a, b}]) -2 = r{G[S U {b}]) = r{G[S]). ■ 



If b — bQ,bi, ...,bk are unlooped and pendant on a then bo,bi, ...,bk are 
identical twins, so we can combine them into a single re- weighted vertex b using 
Theorem [12] and then remove b using Theorem 1131 We refer to the removal of 
any number of unlooped vertices pendant on the same vertex as an unlooped 
pendant vertex reduction. Theorem [13] extends Proposition 4.12 of |2Qj and 
Section 3.2 of Jj, where collections of pendant vertices are called combs. 

A recursive calculation that depends solely on Theorem [7] and Corollary [5] 
is represented by a mixed binary-ternary computation tree; if Corollary [9] is 
always used in place of part (b) of Theorem [7] then the tree will be binary. 
Using identical twin and unlooped pendant vertex reductions makes the com- 
putation tree smaller, because each time one of these reductions is used, we 
avoid splitting the resulting portion of the calculation into branches. Similarly, 
a Tutte polynomial computation that incorporates series-parallel reductions and 
deletion-contraction operations will generally result in a smaller formula than 
a computation that involves only deletion-contraction operations can provide. 
This latter observation was made precise in ,28j: computing the Tutte polyno- 
mial of a matroid M using series-parallel reductions and deletion-contraction 
operations will result in an expression with at least (3{M) terms, where (3{M) is 
Crapo's (3 invariant T8j; moreover if P{M) > the lower bound is attainable. 
(This result is merely the extension to the Tutte polynomial of the important 
theory of reliability domination; see [9l[l5] for expositions.) For interlace poly- 
nomials, analogous lower bounds are derived from the unweighted vertex-nullity 
polynomial q]y{G^). The coefficient of y in q^iG^) is denoted j{G) as in pO] . 
and the evaluation qjv(G")(0) is denoted £{G). 

Theorem 14 If part (c) of Theorem^is applied only to loopless graphs then a 
computation of q{G) using Corollary\^and Theorems^ \lS\ and \13\ is represented 
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by a computation tree with no fewer than ^e(G) leaves. If G is a simple graph 
then the computation tree has no fewer than ^^{G) leaves. 

Proof. We quickly review some elementary properties of qN{G") from [31 
131 [20]. This polynomial is described recursively as follows. If G has a loop at 
a then q^iG^) = qN{{G^Y — o) + qwiG"' — a), if a and b are loopless neighbors 
then gAr(G") = gAr(G" ~ a) + gAr((G")°'' - 6), and if En is the edgeless n- 
vertex graph then gAr (£",") = y". The latter includes the empty graph Eq with 
Qn{Eq) = 1. It follows by induction on the number of vertices that no graph 
has any negative coefficient in (7Ar(G"), and hence every graph has 7(G) > 
and e(G) > 0. Moreover, every nonempty simple graph has e{G) = and every 
disconnected simple graph has 7(G) = 0. 

Theorem [Til is certainly true for as ^j{En), ^e{En) < 1. 

Proceeding inductively, observe that if Corollary [HI or part (a) or (b) of Theo- 
rem [71 is applied to a graph H represented by a certain node of the computation 
tree, and Hi and H2 are the graphs represented by the (first two) resulting 
branch nodes, then = gAr(77{') + Qjv(iJ|); Certainly then ^e{H) = 

^e{Hi) + ie(i?2)- (By the way, we call vertices of the computation tree nodes 
in order to distinguish them from vertices of G.) If Theorem IT^ or [T51 is applied 
to remove a vertex b from a graph H then ^qN{H^){0) — ^qN{H'" — 6)(0), 
because the formulas of Theorems 1121 and [T3l vield a' (a) — (3' [a) = 1 when a (a) 
= /3(a) = a{b) = /3(6) = 1, x = 2 and y = 0. 

It remains to consider the special case involving simple graphs. We actually 
prove a slightly different result, namely: if G is simple then the portion of the 
computation tree involving only nodes corresponding to connected graphs has 
no fewer than ^j{G) leaves. If G is disconnected then 7(G) = so G satisfies 
the result trivially. If G is a connected, simple graph with n < 2 then G satisfies 
the result because 57(G) < 1. Proceeding inductively, observe that if Corollary 
1^1 or Theorem [3 (b) is applied to a connected graph H represented by a certain 
node of the computation tree, and Hi and H2 are the graphs represented by 
the (first two) resulting branch nodes, then (7jv(i/") = q^iHi) + qjv(i?2 )■ It 
follows that l-yiH) ^ \i{Hi) + \'^{H2). If Theorem [H or Theorem [11 is used 
to remove a vertex b from a connected, simple graph H with 3 or more vertices 
then 7(i/) = 7(H - 6), by Corollary 4.17 of [201. ■ 

Theorem [TJl is of limited value because the computations discussed are not 
optimal. The restriction that part (c) of Theorem [71 is only applied to loopless 
graphs is an obvious inefficiency. In addition, if some (combinations of) weights 
are then it would be natural to simply ignore the corresponding parts of the 
computation. There are also other useful twin reductions that do not fall under 
Theorem [HI 

Definition 15 Two vertices a,b of G are fraternal twins if (i) either they are 
looped and nonadjacent or they are unlooped and adjacent, and (ii) they have 
the same neighbors outside {a, 6}. 
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Unlooped fraternal twins are called true twins in [50], but we prefer the 
present terminology because the rows and columns of the adjacency matrix 
corresponding to fraternal twins are not quite the same. Here is an extension 
of Proposition 4.15 of 20] to weighted graphs. 

Theorem 16 Suppose a and b are fraternal twins in G. Let G' he the graph 
obtained from G — b by changing the weights of a: a'{a) — a{a)(3{b) + (3{a)a{b) 
and I3'{a) = f3{a)P{b) + a{a)a{b){x - 1)2. Then q{G) = q{G'). 

Proof. If a e 5 C V{G') then r{G'[S]) = r{G[S]) = r{G[{S U {b}) - {a}]), 
because the adjacency matrices are the same. If a ^ S* C V{G') then r{G'[S]) 
= r{G[S]) =r{G[SU{a,b}])-2. m 

Theorem 1161 has the following inductive generalization. Suppose k > 1 and 
a = bo,bi, ...,bk are fraternal twins in G. If G' is the graph obtained from 
G — fei — ... — 6fe by changing the weights of a to 

«'(«) = E ( n "(^*))( n ^(^M^ - 1)'^'^' 

SC{bo,...,bk} b,eS hj^S 
\S\ odd 

and 13' {a) = ^ ( 11 "(^^))( 11 " 1)'^', 

SC{6o,...,f)fc} fc>6S bj^S 
\S\ even 

then q{G) — q{G'). We refer to the process of combining several fraternal twins 
into a single vertex as a fraternal twin reduction no matter how many fraternal 
twins are combined. 

In general, fraternal twin reductions are just as useful as identical twin re- 
ductions. However, as noted in Proposition 38 of [4] and Corollary 4.16 of [20] 
they have the effect of multiplying the unweighted vertex-nullity polynomial 
by powers of 2. (Observe that if a{a) = /3(a) = a{b) = f3{b) = 1 and x = 2 
then Theorem [TCI gives a' (a) = /3'(a) — 2.) Consequently the lower bounds of 
Theorem [TH arc not valid for computations that utilize Theorem [T51 along with 
Corollary [9] and Theorems H [12] and [H 

If we are given a reduction of a graph G to disconnected vertices using twin 
reductions and unlooped pendant vertex reductions, then Theorems 1121 [T51 and 
1161 describe q{G) in linear time - simply update the vertex weights at each 
step, and at the end refer to part (c) of Theorem [7] If we are not given such 
a reduction then determining whether or not any such reduction exists, and 
finding one if possible, can be accomplished in polynomial time: as in Corollary 
5.3 of [5D], simply search V{G) repeatedly for unlooped degree-one vertices and 
pairs of vertices a, 5 with the same neighbors outside {a, 6}. Hence if G has 
such a reduction then in polynomial time, Theorems I12i 1131 and [16] provide a 
description of q{G) that completely avoids the branching formulas of Corollary 
[n]and parts (a) and (b) of Theorem[7] This observation extends Theorem 6.4 of 
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[20] from gjv to q and from simple graphs that can be analyzed without fraternal 
twin reductions to looped graphs that can be analyzed using all three types of 
reductions: 

Theorem 17 If a graph G can he reduced to a collection oj disconnected vertices 
using unlooped pendant vertex reductions and the two types of twin reductions 
then Theorems \12l \T3\ and \J6\ provide a polynomial-time description of q(G). 

The theorem refers to a description rather than a computation because we 
have ignored the cost of arithmetic operations in the ring R. In case R = Q 
arithmetic operations have low cost, and describing a weighted "polynomial" 
q{G) is the same as computing an evaluation of the unweighted polynomial q{G). 
The full unweighted polynomial may be recovered from several evaluations by 
interpolation, so the theorem provides a polynomial-time computation of the un- 
weighted polynomial. In more complicated rings like y, ai, .., an, Pi, /?„] 
arithmetic operations may be so expensive as to prohibit polynomial-time com- 
putation of entire weighted polynomials. A thorough discussion of these matters 
is given by Courcelle [TB]. 

4 Composition 

In this section we reformulate and extend some results of Arratia, BoUobas 
and Sorkin [4J regarding substituted graphs, using the following version of a 
construction introduced by Cunningham |19j . 

Definition 18 A vertex-weighted graph G is the composition of vertex-weighted 
graphs H and K , G ~ H * K , if the following conditions hold. 

(a) V(H) n V{K) consists of a single unlooped vertex a. 

(h) The vertex a is unweighted in both H and K, i.e., a{a) = 1 = P{a) in 
H and K . 

(c) V{G) = {V{H) U V{K)) — a, and the vertices of G inherit their weights 
from H and K . 

(d) E{G) = E{H) U E{K) U {vw\va e E{G) and aw G E{H)}. 

Requiring a to be unlooped and unweighted in both H and K guarantees 
that no information is lost when we remove a in constructing G. 

Definition [TSl includes several other familiar notions. If a is isolated in H or 
K then H * K is the disjoint union ot H — a and K — a. If a is adjacent to every 
other vertex of H and K then H * K is also denoted {H — a) + {K — a). This 
is traditionally called a "join" but that term has recently been used for general 
compositions (12,. If a is adjacent to every other vertex of H then H * K \s the 
graph obtained by substituting H — a for a in K. If a is adjacent to every other 
vertex of H and H — a is edgeless or complete then the (un) looped vertices of 
H — a are twins in H * K. 

The following observation will be useful. 
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Figure 1: Composition of graphs. 



Proposition 19 Suppose two graphs Ti and T2 are identical except for the 
weights of a single vertex a, and let T be the graph that is identical to both Ti 
and T2 except for ar{a) — (a) + ctr2 ('^) '^'^'^ /^r (a) = f^Ti (a) + /3r2 (cf) • Then 
q{Ti)+q{T2)=q{T). 

Let a be an unweighted vertex of a simple vertex-weighted graph H, and let 
Hi and H2 be the full subgraphs of H ~ a induced by the neighbors (resp. non- 
neighbors) of a. Any weighted interlace polynomial q(H * K) may be analyzed 
in the following way. 

Step 1. Eliminate all edges between vertices oiH—a using pivots and weight- 
changes as in Corollary [HI (In general there will be many different sequences of 
pivots that may be used; this lack of uniqueness is not important in the analysis.) 
The assignments of individual vertices of H to Hi and H2 may change during 
this process, and vertex weights may also change; but these reassignments and 
weight changes will be the same for all graphs K. As in the proof of Proposition 
39 of [3] these pivots will not affect the internal structure of K, because no 
two vertices of K have distinct, nonempty sets of neighbors in H . Here the 
phrase "internal structure" refers to vertex weights, the positions of loops, and 
adjacencies, including adjacencies between a and other vertices of K. 

Step 2. When Step 1 is complete, q{H * K) is expressed as a sum in which 
each summand is the product of an initial multiplying factor and the weighted 
interlace polynomial q{H' * K) of a graph in which every edge is incident on 
a vertex oi K — a. If w G V{H2) then v is isolated, so the only effect of v is 
to multiply that summand by q{{v}) = a{v){y — 1) + P{v). This same effect 
is realized by removing v and incorporating a{v){y — 1) -I- /3(u) into the initial 
multiplying factor, so we may assume that every summand has H2 = 0. In a 
summand with |y(iJ{)| > 1 the vertices of H'l are all identical twins, and may 
be consolidated into a single vertex a using Theorem 1121 In a summand with 
|F(if()| = a single vertex a G V{H'i) may be introduced with a{a) = and 
/3(a) = 1; this will not affect the value of the corresponding summand. As in 
Step 1, these manipulations are the same for all K. 
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Step 3. The weighted interlace polynomial q{H * K) is now expressed as 
a sum in which each summand is the product of an initial multiplying factor 
and a weighted interlace polynomial q{K') in which K' differs from K only in 
the weights of a. In each summand we multiply the a and fi weights of a by 
the initial multiplying factor. This has the effect of multiplying q{K') by that 
factor, so the summand is now simply q{K'). Proposition (TH] tells us that the 
sum is equal to a single weighted interlace polynomial q{K'), where a{a) and 
/3(a) are obtained by adding together the a and f3 weights of a in the various 
summands. 

We deduce that there are weights a{a) and (3{a) that depend only on H and 
a, and have the following property: In every instance of Definition [T51 involving 
H, the interlace polynomial q{H * K) equals q{K'), where K' is obtained from 
K by using a{a) and /3(a) as weights for a. Finding explicit formulas for these 
weights is not difficult. 

Theorem 20 Let H be a vertex-weighted simple graph with an unweighted ver- 
tex a. Then every composition H * K has q{H * K) ~ q{K'), where K' is 
obtained from K by using the following weights for a. 



Proof. With V{K) = {a}, we have q{H - a) = q{H * K) = q{K') = 
/3(a) + {y — l)Q!(a). With K consisting of two adjacent, unlooped, unweighted 
vertices a and v we have q{H) = q{H*K) = q{K') = ((x— 1)^+?/— l)a(a)+y/3(a). 
The stated formulas for a{a) and /3(a) follow. ■ 

In case (a; — 1)^ — (y — 1)^ might be a zero divisor in the ring R, one can avoid 
any difficulty with the formulas of Theorem [501 by first evaluating them in the 
polynomial ring Z[a;, y, ai, .., q;„, /3i, /3„], and then evaluating the resulting 
division-free formulas in R. 

Theorem [50] concerns compositions H * K in which H is simple. If H has 
looped vertices, a similar analysis requires two new steps. 

Step 0. Begin by removing all loops m H — a using local complementation 
as in part (a) of Theorem [T] The result is a description of q{H * K) as a sum 
in which each summand is the product of an initial factor and an interlace 
polynomial q{H' * K) or q(H' * if°), with no loops in H' . 

Step 4. After applying steps 1 - 3 to each of these summands, collect terms 
to obtain a formula 



a{a) 



q{H) - yq{H - a) 
(x-l)2-(y-l)2 



/3(a) 



{{x^l)^ + y-l)q{H-a)-{y 
{x-lf-{y-lf 



mm 



q{H*K) = q{K')+q{{K-)'). 
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In order to distinguish the two terms on the right-hand side we denote by Cc the 
copy of a in K"^. It might seem that we now have to determine four unknowns, 
namely the vertex weights a{a), a{ac), /3(a) and /3{ac) It turns out though that 
these unknowns are not independent. There is an obvious isomorphism between 
{K'Y and JC", and consequently 

q{{Kr)-P{a)q{{Kr-a) 

may be obtained from 

q{{K'^)') - P{aM{K'')' - ac) 
simply by replacing a{ac) with a{a). Theorem [5] tells us that 

q{K') - P{a)q{K' - a) ^ qiiK'T) - P{a)q{{Kr - a), 

so 

q{K') - l3{a)q{K' - a) 

may be obtained from 

q{{K'')') - P{aM{K'')' - a,) 

by replacing a{ac) with a(a). That is, the coefficient of a{a) in q{K') is precisely 
the same as the coefficient of a{ac) in q{{K°')'). It follows that the sum q{K') 
+ q{{K°'y) is unchanged if we replace a{a) by a{a) + a{ac) and replace a{ac) 
byO. 

Theorem 21 Let H he a vertex-weighted graph with an unweighted, unlooped 
vertex a. Then H and a determine weights a{a), /3(a) and f3{ac) such that every 
composition H*K has q{H * K) — q{K') + q^^K"')'), where K' is obtained from 
K by using a{a) and [3{a) as weights for a and (if")' is obtained from K"' by 
using a{ac) — and f3{ac) as weights for a^, the copy of a in [K"")' . 

Formulas for the three weights mentioned in Theorem [^l] may be derived 
from three instances of the theorem. We use H — a and H as in Theorem 
[201 and also H^, the graph obtained from H by attaching a loop at a. These 
correspond respectively to compositions of H with graphs Ki, if 2, such that 
V{Ki) = {a}; V{K2) = {a,v} with v an unweighted, unlooped neighbor of a; 
and ^(ifa) = {a,i} with £ an unweighted, looped neighbor of a. Definition [1] 
gives the following values. 

qiH - a) = q{K[) + q{{Kl)') ^ {y - l)a(a) + /3(a) + /3(ae) 



q{H) = q{K'^) + q{{K^)') - {{x - if ^ y ~ l)«(a) + yP{a) + xP{a,) 



q{H') = q{K'^) + q{{K^)') = {{x - 1)^ + y - l)a(a) + xP{a) + y/3(a,) 
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We deduce these formulas. 

(x + y)q{H -a)- q{H) - q{H') 



a{a) 



(a; + 2/)(y-l)-2((.T-l)2 + y-l) 
q{H) - q{H') 



p{a) - /3(a,) 



y- X 



2((x - 1)^ + y - 1) - (y - l)[x + y) 

Separate formulas for (3{a) and (3(ac) are derived in the obvious ways by adding 
and subtracting the last two. As before, possible problems with denominators 
may be avoided by first evaluating the formulas in Z[a;, y, ai, .., q;„, /3„]. 

We close this section with different formulas for the weights a(a), (3{a) and 
/3(ac) that appear in Theorems [20] and [21] Suppose H has an unweighted, 
unlooped vertex a. Let iV(a) denote the open neighborhood of a, i.e., the 
set containing the vertices v ^ a E V{H) that neighbor a. Given a subset 
S C V{H — a) let p — ps,a (resp. k = Ks,a) be the row (resp. column) vector 
with entries indexed by {i : Vi G S} whose i"* entry is 1 or according to 
whether Vi € N{A) or Vi N{a). Also let AI — Ms be the adjacency matrix of 
H[S]. Note that r{M) < r (M k) < r{M) + 1 and 

r(Af) <r ^ j ' p 1 j 

because adjoining a single row or column to a matrix raises the rank by or 1. 
Definition 22 The type of S with respect to a is defined as follows. 

S is of type 1 if r{AI) = ^f^^ o)^''(^^ l)^"^' 



S is of type 2 if r{M) 



S is of type 3 if r{M) = r 



AI k\_ ( M k 
P j p 1 



AI K \ _ ( M K 

p 1 j P 



Lemma 23 Every S C V{H — a) is of one of these types. Moreover, if a has 
no looped neighbor then there is no S V{H — a) of type 3. 

Proof. The fact that every 5 C — a) is of type 1, 2 or 3 appears in 

Lemma 2 of [5]- 
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Suppose S is of type 3, so 



"(M) + l = r 



M K 

9 



The row vector {p O) cannot be a sum of rows of (A/ k) . For if it were then 
p would be the corresponding sum of rows of M , and by symmetry k would be 
the corresponding sum of columns of M. Consequently it would follow that 



-(M) = r{ M 



K ) = r 



M K 

P 



On the other hand, p must be the sum of the rows of M corresponding to the 
elements of some subset T ^ S, for if it were not then it would follow that 

r(Af) + l=r(^^^)=r( Af O = - ( "Jf l)-^- 

Every such T must contain an odd number of neighbors of a, to avoid giving a 
sum of rows of (A/ k) equal to (p O) . 

Choose such a T, and consider the induced subgraph H[T n N{a)\. As the 
sum of the rows of M corresponding to elements of T is p, the sum of the rows 
of the adjacency matrix of H\T n N{a)] is the row vector (1...1). |T n N{a)\ is 
odd, so the adjacency matrix of H\Tr\N{a)\ has an odd number of entries equal 
to 1. The matrix is symmetric, so at least one of these entries must appear on 
the diagonal. That is, at least one vertex of H[T n N{a)] is looped. ■ 

For S C V{H — a) let Af" = Mg be the matrix obtained from M by toggling 
every entry rriij that has Vi,Vj S N{a). Then r (Af — r (Af Kj because 
the first matrix is transformed into the second by adding the last column to every 
column corresponding to a neighbor of a. Similarly, adding the last column of 
the first matrix below to every column corresponding to a neighbor of a tells us 
that 

' M K \ _ f M'' K 

P I 1 



1 + r(Ar ) 



Consequently, Definition [22] may be restated using the relationship between 
r{M) and r(Af"): if S is of type 1 then r(Af ) = r{M''), if S is of type 2 then 
r\M) ^ ri^M") ~ 1, and if S is of type 3 then r{M) = r(Af'^) + 1. 



Proposition 24 Fori e {1,2,3} let 

S(lV(H-a) seS v^S 
of type i 



n{(H^a)[S\) 
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Then the following equations hold. 

q{H -a)=qi{H-a) + q2{H - a) + qsiH - a) 



q{H) ^yq^{H -a)+[l + l^Jll \ q^lH - a) + xq^iH - a) 



q{H') =xq^{H -a)+[l + ) q2{H - a) + vq:,{H - a) 



Also, 



q{H'' -a) = q,{H-a)+ 92(i/ - a) + (f— Q^iH - a). 

Proof. The first equality is obvious. For the second, note that each S C 
V{H — a) gives rise to two subsets of V{H), namely S and S U {a}; these 

correspond to the adjacency matrices M and ( ) ' ^^^^^'^^y' third 

equality is derived by considering the contributions of two adjacency matrices 
for each S C V{H — a), namely M and i ) ■ The last equality follows 



P 1, 

from the discussion preceding the proposition. ■ 

Corollary 25 The weights mentioned in Theorems \2U\ and \21\ are /3(a) = qi{H— 
a), a(a) = q2(,H — a)/(y — 1) and f3{ac) — qz{H — a). In particular, (3(ac) — 
if a has no looped neighbor. 

Proof. The corollary follows from Proposition [21] and the formulas given in 
Theorem BUI and immediately after Theorem [5TJ ■ 

5 A characterization of simple graphs 

In this section we focus our attention on the unweighted vertex-nullity polyno- 
mial. 

Proposition 26 If G is a connected, unweighted graph with at least one looped 
vertex then e{G) = gjv (G)(0) > 1. 

Proof. The proposition is certainly true if G has n < 2 vertices, as all three 
such graphs have e(G) — 2. The argument proceeds by induction on n > 3. 
Recall that e{G) > for every graph G, and let a be a looped vertex of G. 
Then qN{G) = qN{G - a) + qniG" - a), so e{G) > max{£(G - a), e(G" - a)}. 

Suppose a is not a cutpoint of G. If G has some looped vertex other than a 
then the inductive hypothesis implies that e(G — a) > 1. If G has no looped ver- 
tex other than a then every neighbor of a is looped in G" — a. Every component 
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of — a contains at least one neighbor of a, so the inductive hypothesis imphes 
that e{C) > 1 for every component C of G°-—a. Hence e{G°'—a) = He ^(^) 1- 
Suppose now that G has a looped cutpoint a. For each component C of G — a 
and each vertex v G N{a) that lies in some other component of G — a, no edge 
connecting v to an element of N{a)nV{C) appears in G — a. Consequently every 
such edge appears in G° — a, so all of N{a) is contained in a single component 
of G" — a. As G is connected, this implies that G" — a is also connected. Hence 
if any neighbor of a is unlooped in G, G° — a is a connected graph with a looped 
vertex and the inductive hypothesis implies that e(G° — a) > 1. If instead every 
neighbor of a is looped in G, then every component G of G — a has a looped 
vertex and the inductive hypothesis implies e{G — a) — Ylc ^(C*) > 1- ■ 

Corollary 27 An unweighted graph G has e{G) > if and only if G has no 
nonempty simple component. 

Proof. As noted in Remark 20 of [J, every nonempty simple graph has 
e(G) = 0. It follows that every graph with a nonempty simple component also 
has e{G) — 0. On the other hand, Proposition and the fact that e{Eo) = 
1 together imply that every graph with no nonempty simple component has 
e(G) > 0. ■ 

Corollary 28 Let G be an unweighted graph, and let G^ be its complement, i.e., 
the graph with V{G'^) ~ V{G) whose edges (including loops) are precisely the 
edges absent in G. Then the following are equivalent: G is simple, g7v(G+£'") = 
yqN{G') andy\qNiG + E^). 

Proof. Suppose first that G is simple, and let ff = G + i?" with a the 
vertex of i?". Let K be the two-vertex graph with V{K) — {a, £} in which £ is 
an unweighted, looped neighbor of a. Theorem [20l tells us that 

q{H *K)^ q{K') ^ xP{a) + {{x - 1)^ + y - l)a(a). 

On the other hand, the recursive description of q [5 tells us that 

q{H * K) = q{{H *K)-£) + {x- l)q{{H * KY -£) = q{G) + (x - 1)9(G"). 

Recalling that q{G) = q{H — a) = (3{a) + {y — l)a{a)^ we see that 

= q{H*K)- q{H * K) ^ {x - l)/3(a) + {x - lfa{a) - [x - VjqiG") 

and consequently g(G^) — (3{a) + {x — l)a{a). Recalling that g(G + £'}') ~ q{H) 
= {{x — 1)^ + y — l)a(a) + yf3{a), we see that 

gjv(G + E^) = q{G + £;nU=2 = y {a{a) + /3(a)) U=2 = yq{G^)\x=2 = 2/<ZAr(G^). 

Now suppose that G is not simple. Then G has a looped vertex, and so 
does G + E]". As G + E]" is connected, Corollary [27] states that e{G + E^) > 0; 
consequently 2/ f (7Ar(G + -E"). ■ 
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6 Trees 



A combinatorial description of the interlace polynomials of trees and forests is 
given in [2j. In order to motivate the next section we sketch this description 
briefly here, omitting details and proofs. Recall that a tree T is rooted by 
specifying a root vertex r G V{T). Each non-root vertex v G V{T) then has a 
unique parent p(v), a neighbor whose distance from r is less than the distance 
from r to v. The elements of p~^{{p{v)}) are the children of p(v), and the 
children of p{v) other than v itself are siblings of v. An ordered tree is a rooted 
tree given with an order on the set of children of each parent vertex; non-root 
vertices may then have earlier siblings and later siblings. A set of vertices that 
contains no adjacent pair is independent, and a set of vertices dominates a vertex 

V if it contains v or contains some neighbor of v. 

Definition 29 An earlier sibling cover (or es-coverj in an ordered tree T is an 
independent set I that dominates r and has the property that for every non-root 
vertex v G I, every earlier sibling of v is dominated by I. 

Definition 30 For integers s and t the es-number Cs^t{T) is the number of 
s-element es-covers in T whose non-root elements have t different parents. 

If T' is a subtree of T and r e V{T') then we presume that the children of 
each parent vertex v in T' are ordered by restricting the order of the children of 

V in T. With this convention, it is easy to verify that the earlier sibling covers 
in large trees arise from earlier sibling covers in subtrees. 

Lemma 31 Let T be an ordered tree with a leaf £ such that p{£) ^ r ^ I, all 
the siblings of i are leaves, and £ has no later siblings. Then 



{es-covers I in T with £ £ 1} 
= {unions p~^{{p{£)}) U / with I an es-cover in T — p{£) — P^^{{p{£)})}- 

Lemma |31] is the key to an inductive proof that the terms in the definition 
of q{T) can be collected into sub-totals corresponding to earlier sibling covers. 

Definition 32 Let T be an ordered tree with vertex weights, and let I be an 
es-cover in T. Let 



{es-covers I in T with £ ^ 1} ^ {es-covers in T — £} 



and 



Ir 
h 



{v e I 
{vel 

{v^I 



either v = r or I contains a later sibling of v}, 
V ^ r and I contains no later sibling of v}, and 
I contains a child of v} 
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For each vertex v define the /-weight wi{v) as follows: 

( Piv) + a{v)iy ~ 1) ifveir 

^^^^^ ^ \a{v) ■ {{y - + - ^?o^{p{y))) ^fv^h 

I 1 ifvel'^ 

^f3{v) if V ^ I and V ^ I'^. 

The product 

n ^'f(^) 

vGV(T) 

is the total weight of I in T , denoted wt (I) ■ 

Theorem 33 // T is an ordered tree with vertex weights then 

q{T) = ^^W- 
/ an es — cover 

Proof. If T has no more than one vertex of degree > 2, it is not difhcuh to 
verify the theorem directfy. If T has more than one vertex of degree > 2 then 
the theorem follows inductively from Theorem [7] and Corollary [5] using Lemma 
[511 as detailed in [2]. ■ 

Setting a = 1 and /3 = 1 we deduce that the unweighted interlace polynomial 
of a tree is determined by a simple formula involving earlier sibling covers. 

Corollary 34 If T is a tree then the unweighted interlace polynomial of T is 

s,t 

As the interlace polynomials are multiplicative on disjoint unions, these re- 
sults extend directly to disconnected forests. 



7 Algorithmic activities 

The Tutte polynomial is a very useful invariant of graphs and matroids, which 
incorporates a great deal of information and can be defined in several different 
ways; see [10], [29] and [30] for more detailed discussions than we provide here. 
One of its definitions is given in Definitions [35| and [36l 

Definition 35 Suppose G is an unweighted graph with E{G) = {ei,...,em}, 

and T is a maximal spanning forest of G. An element Ci ^ E{T) is externally 
active with respect to T if i is the least index of an element of the unique circuit 
contained in E{T) U {e,}. An element Ci G E(T) is internally active with 
respect to T if i is the least index of an element of the unique cutset contained 
in {E{G) — E{T))U{ei} . The numbers of edges that are externally and internally 
active with respect to T are denoted e{T) and i{T), respectively. 
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Definition 36 If G is an unweighted graph with E{G) = {ei, ...,e„i} then the 
Tutte polynomial of G is 

t{G)=Y,x^^^)y<^\ 

T 

Theorem [33] and Corollary [M] bear a strong resemblance to this formula, 
with earlier sibling covers in rooted trees replacing maximal spanning forests in 
arbitrary graphs. We do not know whether there is a combinatorial analogue 
of Definition [36] that describes the interlace polynomials of an arbitrary graph, 
but there is an algorithmic analogue. Before presenting it, we recall another 
definition of the Tutte polynomial. 

Definition 37 IfG is an unweighted graph then the Tutte polynomial t{G) is an 
element of the polynomial ring l^x, y\ determined recursively by these properties. 

(a) If e G E{G) is neither a loop nor an isthmus then t{G) = t{G — e) + 
t{G/e). 

(b) IfXe E{G) is a loop then t{G) = yt{G - A). 

(c) If(3£ E{G) IS an isthmus of G then t{G) = xt{G/f3). 

(d) For any positive integer n, t{En) ~ 1. 

To recursively calculate t{G) one simply chooses an arbitrary edge of G, and 
applies the appropriate part of Definition [37| this process is repeated as many 
times as necessary. Such a computation is represented by a computation tree 
in which a node that represents an instance of part (a) has two children and a 
node that represents an instance of part (b) or (c) has only one. 

Proposition 38 Let G be a graph with E(G) = {ei, ...,em}, and consider the 
recursive implementation of Definition \37\ in which Cm is removed first, then 
Cm-i is removed in all branches, then 6^-2 is removed in all branches, and 
so on. The leaves of the computation tree representing this implementation 
correspond to the maximal spanning forests of G, with the leaf corresponding 
to T resulting from the portion of the computation in which edges of T are 
contracted and elements of E{G) — E{T) are deleted. A node of this portion of 
the computation tree represents the removal of an active edge if and only if it 
has precisely one child. 

Proof. The proposition is implicit in the fact that Definition [37l and Defini- 
tion[3n]both yield t{G), so it appears implicitly in just about every presentation 
of the Tutte polynomial. See Theorem IX. 65 of [29 or Theorem X.IO of [TD], 
for instance. Explicit discussions of the connection between activities and com- 
putation are less common in the literature, though there are some j6l [21] [22]. 

The proof is a direct induction on |_B(G)|. If e„i is a loop then the maximal 
spanning forests of G and G — Cm coincide, and the computation tree for t{G) 
is obtained from the computation tree for t{G — &,„) by attaching a new root 
node of degree 1, representing the removal of e™. If e™ is an isthmus then 
the maximal spanning forests of G and G/em correspond, and the computation 
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tree for t{G) is obtained from the computation tree for t{G/em) by attaching 
a new root node of degree 1. Otherwise, the maximal spanning forests of G 
that contain Cm correspond to the maximal spanning forests of G/em, and the 
maximal spanning forests of G that do not contain Cm are the maximal spanning 
forests of G — Cm- The computation tree for t{G) consists of the root and two 
disjoint subtrees that are the computation trees for t{G/em) and t{G — em). ■ 

Definition [35l defines active edges using the structure of G, and this leads to 
Definition[36ls description of t{G) as a generating function for maximal spanning 
forests. Proposition [38] shows that we may also see activity from an algorithmic 
viewpoint: the external activity of a particular e ^ E{T) is revealed in the 
fact that one step of a calculation of t{G) involves removing e using part (b) 
of Definition [571 rather than part (a). This distinction affects the result of the 
computation, so it would be important even if activity could not be conveniently 
described using the structure of G, or did not contribute to a convenient closed 
form for t{G). 

Here is an analogue of Definition [37] for the weighted interlace polynomial. 

Definition 39 If G is a weighted graph then q{G) is determined recursively by 
the following properties. 

(a) If a is a looped vertex then 

q{G) = P{a)q{G -a)+ a{a){x - l)q{G'' - a). 

(b) If a and b are loopless neighbors in G then 

q{G) - p{a)q{G - a) + a{a)q{{G-' - 6)'), 

where [G""^ — b)' is obtained from G"''^ ~b by changing the weights of a to a' [a) = 
(3{b) and P'{a) = a{b){x ~ 1)2. 

(c) If a is isolated and looped then q{G) — (a(a)(x — 1) + j3[a))q{G — a). 

(d) If a is isolated and unlooped then q{G) — {a{a){y — 1) + P{a))q{G — a). 

(e) The empty graph has q{%) = 1. 

The preceding discussion of activities and the Tuttc polynomial suggests the 
following. 

Definition 40 A node of a computation tree representing a recursive implemen- 
tation of Definition \39\ is active if it has precisely one child, i.e., if it represents 
an application of part (c) or part (d). 

An "activities formula" for q{G) arises directly from a computation tree 
representing an implementation of Definition[39] The formula has one summand 
for each leaf of the computation tree (each call to part (e) of Definition [5^ , 
representing the product of the coefficients contributed by the nodes in the 
portion of the computation tree that gives rise to that leaf. Active and non- 
active nodes contribute different coefficients. 
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If G is a rooted tree then the formulas of Theorem [551 and its corollaries are 
activities formulas. 

Proposition 41 Let T be a rooted tree with root r, and consider a recursive 
implementation of Definition \39\ structured as follows. If possible, apply part (d) 
of Definition \39[ if not and there is a parent vertex other than r then apply part 
(b) with a leaf of the type denoted I in Lemma\3^ as a; otherwise apply part (b) 
with the last child of r as a. The leaves of the computation tree representing this 
implementation correspond to the earlier sibling covers of T , with the es- cover 
I corresponding to a given leaf constructed from the portion of the computation 
that gives rise to that leaf as follows: an occurrence of part (d) of Definition \39\ 
contributes its a to I, and an occurrence of the q{{G°''' — b)') branch of part (b) 
contributes its a to I . 

Proof. If T has no vertex other than r then {r} is the only es-cover in T, 
and the computation consists simply of a single call to part (e) of Definition [39l 
If T contains no parent vertex other than r and a is the last child of r, then 
the es-covers in T include V{T) — {r} and the es-covers of T — a. The first step 
of the computation is an application of part (b) of Definition [551 with b = r. 
The inductive hypothesis applies to T — a, and the branch of the computation 
corresponding to (T"'' — 5)' = (T — r)' consists solely of calls to part (d) because 
every vertex of T — r is isolated; consequently the latter part of the computation 
tree contains only one leaf, corresponding to V{T) — {r}. If T contains a parent 
vertex other than r, then the first step of the computation is an application 
of part (b) of Definition [39l with a = ^ as in Lemma [3T1 and b = p{£). The 
computation tree contains a single node representing this first step and also two 
subtrees, one corresponding to T — a = T — £ and the other corresponding to 
(T"^ — b)' = {T — p{£)y . The proposition follows inductively from Lemma [511 
■ 

We do not know whether or not it is possible to reformulate Definition [40l so 
that it always refers to G instead of a computation tree. Such a reformulation 
would certainly be valuable, as the resulting activities formulas would shed light 
on the combinatorial significance of the interlace polynomials. 

Definition [40l extends directly to computation trees representing implemen- 
tations of other recursions. For instance, if Definition [39l is augmented by in- 
corporating pendant-twin reductions then the resulting computation trees will 
have active nodes representing these reductions, in addition to active nodes 
representing parts (c) and (d) of Definition [551 
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